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Abstract 


A new electrode model involving the fractal structure of activated carbon used as electrode material was proposed for an electric double layer 
capacitor (EDLC). The fractal structure of activate carbon was simulated by branch pore structure of three sizes of cylindrical pores. Three sizes 
of cylindrical pores were related to macro, meso, and micro pore, since the pore size of activated carbon has wide distribution with a few modes of 
nm order. The impedance spectrum of EDLC describes the locus of blocking electrode in low frequency range, and the curve at an angle of almost 
45° to real axis in high frequency range on the Nyquist plane. The low and high frequency ranges of the impedance spectrum were defined as a 
lumped constant-type and a distributed constant-type, respectively. Computer simulation of electrochemical impedance with the present electrode 
model was carried out to understand the relation between the impedance and the electrode structure. The contributions of five parameters to 
impedance spectrum were discussed, i.e., depth of pore, diameter of pore, specific resistance, the interfacial impedance at electrode/solution, and 
branch number. The specific resistance p influenced on the shape of impedance spectrum in distributed constant-type range significantly. On other 
hand, the interfacial impedance at electrode/solution interface controlled the shape of the impedance spectrum in lumped constant-type range. In 
the course of the curve-fitting to impedance spectrum of EDLC, the separated impedance spectra related to macro, meso, and micro pores were 
obtained, and the roles of these pores on electric capacity were discussed. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


An electric double layer capacitor (EDLC) is an energy con- 
version device that stores electric charge in the electric double 
layer at electrode/electrolyte interface. Activated carbon with 
high surface area is usually used as the electrode material since 
the usable charge stored in EDLC mainly depends on the surface 
area. In addition, the usable charge is also related to the pore size 
and the pore structure. It was reported that the feature of high 
capacitance in EDLC depended on the characteristics of pore of 
nm order (micro pore) [1-5]. If the ion radius of electrolyte is 
large comparing with the pore radius, the ions cannot penetrate 
completely into the micro pore and all surface area cannot work 
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to store the charge [1,6—9]. Therefore, it can be considered that 
the pore structure of activated carbon is one of the important fac- 
tors for capacitance of EDLC. For example, the roles of macro, 
meso, and micro pores are different to store the charge. 
Electrochemical impedance spectroscopy (EIS) allows the 
information concerning electrochemical and electric properties 
of electrode and electrolyte. Several time constants in electrode 
reactions can be discriminated, and the electrode structure can 
be analyzed by EIS. Therefore, many researchers have used EIS 
to analyze structure of porous electrode [9-18]. Some groups 
adopted transmission line model (TLM) to simulate the fre- 
quency response in a pore of EDLC. For example, Qu and Shi 
[6] carried out the curve-fitting to impedance spectrum of EDLC 
by the assumptions of cylindrical pore with distributed capaci- 
tance and resistance in TLM. They [6] discussed the relationship 
between time constant of charge/discharge and cell voltage. 
Song et al. [19,20] simulated impedance spectrum of EDLC 
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Nomenclature 


Parameter list of single pore 


cross-section area of a pore (cm?) 

capacitance per unit area (F cm~?) 

depth per segment (cm), dX = X/N 

diameter of a pore (cm) 

segment number 

number of partition (piece) (n = 1-N) 

reaction resistance per unit area (Q cm?) 
solution resistance ({2) 

solution resistance per one segment (Q), Rsolseg = 
Ryo dX 

solution resistance per unit length (Qcm7!), 
Roo = p/A 


depth of a pore (cm) 
impedance of a pore (Q) 


impedance per one segment (Q), Zseg = 
Z* /dX (b) 
impedance per unit length (Qcm), Z* = 
Z** /D (a) 


impedance per unit area (Q cm?) 


Greek letters 


specific resistance (Q cm) 


Parameter list of a pore with repeated structure 


cross-section area of a pore per segment of seg- 
ment number m (cm?) 

branch number 

capacitance per unit area (F cm~?) 

depth of a pore per one segment of segment num- 
ber m (cm), dX, = Xm/N 

diameter per one segment of segment number m 
(cm) 

segment number (n = 1-N) 

reaction resistance per unit area (Q cm?) 
solution resistance of segment number m (2) 
solution resistance per unit length of segment 


number m (Qcm7!), Raia = 0/Am 

solution resistance per one segment of segment 
number m (Q), Rsol,m = ReoimoXm 

depth of a pore of segment number m (cm) 
impedance of segment number m (2) 

impedance per unit length of segment number m 
(Qem), ZF, = Z**/Dm 

impedance per one segment of segment number 
m (&), Zseg,m = Za [dX m 

total impedance of pore electrode (2) 
impedance per unit area (Q cm?) 


Greek letters 


specific resistance (Q2 cm) 


macro pore 
p meso pore 


micro pore 


Fig. 1. Conceptual scheme of pore of activated carbon. 


by assuming cylindrical pore electrode, and proposed TLM with 
pore size distribution (TLM-PSD). The simulated results of elec- 
trochemical impedance with these models [6,19-20] are in good 
agreement with experimental ones geometrically. However, the 
pores of activated carbon have complex branch structure as 
described in Fig. 1, and it is considered that impedance anal- 
ysis model should involve the branch pore structure. 

In the present paper, we propose a new conceptual model 
of EDLC with branch pore structure. The single pore electrode 
is represented by a cylindrical pore. The branch pore structure 
of activated carbon was simulated by the repetition with three 
sizes of cylindrical pores in this model. These pores are related 
to macro, meso, and micro pore in descending order. A great 
number of branch pores are located on the electrode, and the 
electrochemical impedance is assessed by sum of the branch 
pores. Computer simulation of electrochemical impedance with 
the present electrode model was carried out to characterize the 
impedance spectrum of EDLC. Furthermore, the contributions 
of macro, meso, and micro pores to electrochemical impedance 
of EDLC were studied by the curve fitting to impedance spec- 
trum. 


2. Theory 


2.1. Lumped constant-type and distributed constant-type 
equivalent circuits 


The schemes of blocking electrochemical impedance that 
excludes Faradaic process are shown in Fig. 2 for (a) lumped 
constant-type equivalent circuit and (b) distributed constant-type 
equivalent circuit. In the present paper, the lumped constant- 
type equivalent circuit can be expressed as a series of solution 


(a) (b) 


lumped constant-type equivalent circuit (a) 


(MHE ) 


ia distributed constant-type equivalent circuit (b) 
(see Fig. 3.) 


Z"/Q. 


distribution 


Reo 
ZIQ 


Fig. 2. Schemes of electrochemical impedance of EDLC. (a) Lumped constants 
equivalent circuit and (b) distributed constants equivalent circuit. 
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Fig. 3. Distributed constant equivalent circuit of a single pore. 


resistance and capacitance. The distributed constant-type equiv- 
alent circuit is composed of solution resistance and impedance 
per one segment as depicted in Fig. 3. The impedance of 
equivalent circuit (a) describes locus in parallel with the imag- 
inary axis on complex plane. The real part is solution resis- 
tance Reo, and the imaginary part converges to zero with high 
frequency. 

The impedance of equivalent circuit (b) has almost similar 
shape to that of the equivalent circuit (a) in low frequency range. 
And it describes the locus of 45° to real axis in high frequency 
range. The reason is that the current distribution on the electrode 
is non-uniform in high frequency range. 


2.2. Theoretical formula of impedance with cylindrical 
electrode by TLM 


Levie [21] solved theoretical formula of TLM impedance for 
inner wall of cylindrical electrode as follows: 


[Re X 
R*,Z* coth sa (1) 


Zz" = i (2) 
mD 
p 
a= 8) 


where X is the depth of a pore, D the diameter of a pore, Z“* the 
impedance per unit area, Z the impedance per unit length, Z the 
impedance of a pore, Ry the solution resistance per unit length, 
p the specific resistance, and A is the cross-section area of a 
pore. Z“* including Faradaic process such as redox capacitor is 
represented by formula (4), and Z** excluding Faradaic process 
(blocking electrode) is represented by formula (5). 


R*t* 
1 +joR* C 
1 
Te = joc (5) 


jis the imaginary number unit, w the angular frequency (=27f, f 
is frequency), R¢* the reaction resistance per unit area, and C3; 
is the capacitance per unit area. 


2.3. Impedance simulation of single pore by using matrix 


[22] 


In the present paper, impedance spectrum is basically simu- 
lated by assuming cylindrical pore as depicted in Fig. 4(a). The 
cylindrical pore is divided into N segments in the direction of the 
depth, and each segment is defined as n= 1, 2, ..., N from top 
to bottom of the pore. Impedance per unit area Z" and solution 
resistance per unit length R¥ were converted to impedance per 
one segment Zseg and solution resistance per one segment Rgolseg 
by formulas (6)-(8), 


Zs 

Zseg = dx (6) 
Vii 

ZY = =D (7) 

Rgolseg = R% dX (8) 


where dX is depth of a pore per one segment (dX =X/N). An 
equation related to current flow in the electrolyte and at the elec- 
trode surface by Kirchhoff’s voltage law was derived in order 
to express the impedance of pore electrode by TLM. There are 
the mesh current method and the branch current method for the 
calculation method by Kirchhoff’s voltage law. The mesh cur- 
rent method was employed in this analysis because it is easy 
to express the equations by a matrix. Fig. 3 shows the rela- 
tion between TLM and the current loops, J; is total current, J, 
(k=2-n) is current of loop k. The equations of the potential drops 
for the current loops in TLM equivalent circuit are summarized 


(a) 


Fig. 4. (a and b) Conceptual model for single pore electrode. 
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as follows: 


N-1 N-1 
E= [E ras + Zas) h+ X Rsolseg In 


n=1 n=2 
+ (Rsolseg F Zseg) IN 

0 = Reotsegli + { Rsolseg + 2Zseg } I2 — Zseg l3 

0= Rsolseg l1 B Zsegln—-1 + { Rsolseg + 2Zseg } In (9) 
m Zsegln+1 

(n = 2,3,...,N—1) 

0 = { Reolseg + 2Zseg} T1 — Zsegln—1 
+ {Rsolseg + 2Zseg } In 


where Fis the potential modulation of the electrode. These equa- 
tions are expressed by a matrix M as follows: 


E h 
0 h 
0 i 
=M| (10) 
0 In-1 
0 In 
where, 
N-1 
5 Rsolseg F Zseg Rsolseg Rsolseg 
N=1 
Rsolseg Rsolseg +2Zso è —Zso 0 
M = | Rsolseg — Zsol 
0 
Rsolseg : 
Rsolseg + Zsol 0 nae 0 


The currents /;—Jy can be assessed by using inverse matrix M~! 
of Eq. (11). 


h E 

h 0 
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=m | a2) 

In-1 0 

In 0 


Impedance Z can be calculated from Eq. (13). 


E 


Z = — 
h 


(13) 

In order to confirm the suitability of the algorithm and the 
calculation method of the present simulation, the impedance 
calculated by Eqs. (12) and (13) is shown in Fig. 5 with the results 
calculated by theoretical expression (1). Both plots are perfectly 
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Fig. 5. Calculated result of electrochemical impedance of single pore elec- 
trode. (O) Calculation with the present method using matrix, (x) calculation 
according to Eq. (1); (X) 75 nm, (D) 50nm, (C%") 2.0 x 10-4 F/cm?, and (p) 
5.0 x 108 Wem. 


Rsolseg + Zsol 


0 


(11) 
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—Zsol 
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identical, indicating that impedance of distributed constant- 
typed equivalent circuit can be calculated by using the matrix. 
In this calculation, the segment number N was 200. The soft- 
ware for the matrix calculation was Matlab. Though the values 
of Rž p Ré*, and C3 were constant to calculate the impedance 
in Fig. 5, the impedance of pore electrode which has other shape 
except for cylinder can be calculated by the same algorithm. 


2.4. Impedance simulation of pore electrode with repeated 
structure 


The electrode model with repeated structure is proposed to 
express the fractal structure of porous carbon electrode. A single 
pore of activated carbon is simulated by the cylindrical pore, 
and the fractal structure of activated carbon illustrated in Fig. 1 
is expressed by the electrode model with repeated structure as 
depicted in Fig. 6(a). The secondary pore is located at bottom 
of the primary pore, and the third pore is located at the bottom 
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(a) 


(b) Rsol, 2 


(c) Rsol, 1 Rol, 2 


Fig. 6. (a) Conceptual model for pore electrode with repeated structure, (b) 
equivalent circuit for pore electrode with repeated structure, and (c) simplified 
equivalent circuit of pore electrode with repeated structure. 


of the secondary pore subsequently. The equivalent circuit for 
the repeated pore structure is shown in Fig. 6(b), where Zm is 
the impedance of mth pore. The solution resistance Rsol,n of mth 
pore is expressed by Eq. (14), 


N 
Rsol,m = X Rsolseg,m (14) 
n=1 


where Rgolsegym Means solution resistance per one segment of 
mth pore. The total impedance of pore electrode with repeated 
structure Ziotal is expressed by Eq. (15). 


1 1 1 


Zui A Raa iG) aye Rate) 


(15) 


The equivalent circuits of three kinds of pores are shown in 
Fig. 6(c). The impedance of macro pore Zma is identical to that 


of primary pore Z1. 

Zma = Z1 (16) 
The impedance of meso pore Zme is series Of Rgoj,1 and Zp. 
Zme = Z2 + Rsol,1 (17) 


In addition the impedance of micro pore Zmi is series of Rsol,1, 
Rsol,2, and Z3. 


Zmi = Z3 + Rsol,1 + Rsol,2 (18) 


Branch number of meso pore against macro pore is defined as 
bı and that of micro pore against meso pore is defined as b2. In 
case that bı or b2 is larger than unity, Ziotal is expressed by Eq. 
(19). 


1 1 1 
=>+ 
Ztota Zt Root + 1/(b1(A/Z2) + (1 /(Rso,2 + 62Z3)))) 
(19) 

And 

Zma = Z1 (20) 
Z2 

Zme = Rsol,1 api (21) 
by 
Rgo1,2 + (Z3/b2) 

Zmi = Rsol,1 + Le (22) 


by 


For example, the scheme of electrode model of brunch pore 
structure with b=3 and the equivalent circuit are shown in 
Fig. 7(a and b). This model has the structure in which three meso 
pores exist at the bottom of one macro pore and three micro pores 
also exist at the bottom of one meso pore. The equivalent circuits 
of three kinds of pores are shown in Fig. 7(c). The impedance 
of macro pore Zma is the same as Z;. The impedance of meso 
pore Zme is composed of Rsoj,1 and three Z2. Furthermore, the 
impedance of micro pore Zmi is composed of Rgo1,1, three Rsol,2, 
and nine Z3. 


3. Experimental method 
3.1. Electrochemical impedance measurement 


EDLC by HONDA R&D (Ultra capacitor) shown in Fig. 8 
was used for electrochemical impedance measurement. Poten- 
tiostat (Hokuto, HG-501G) and frequency response analyzer 
(Solartron, SI 1254) were equipped to measure the impedance 
spectra. The frequency range of impedance was from 10 mHz to 
1 kHz, and the frequency of imposed signal was scanned loga- 
rithmically at 5 point for each decade from the high frequency 
region. The measurements were carried out by galvanostat mode 
and the amplitude of the sinusoidal input current was controlled 
to make the amplitude of ac potential response smaller than 
10 mV. 
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Fig. 7. (a—c) Conceptual model of pore electrode with branch structure. 


4. Results and discussion 
4.1. Electrochemical impedance measurement of EDLC 


Experimental result of electrochemical impedance of EDLC 
at the bias voltage of 2.0V is shown in Fig. 9. Impedance 
spectrum can be divided into two frequency ranges, namely, 
the distributed constant range above 100 mHz and the lumped 
constant range below 100 mHz. The real part of impedance is 


constant, and the imaginary part increases with the decrease of 
the frequency at lumped constant range, indicating the block- 
ing electrode behavior in the low frequency range. Impedance 
deviates from the locus of blocking electrode at distributed 
constant range. It is considered that the current distribution is 
non-uniform by porous structure of electrode in the high fre- 
quency range. The solution resistance at high frequency limit 
is approximately 0.002 Q, indicating the low contact resistance 
and electrolyte resistance of the present EDLC. 
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Fig. 8. Photograph of scheme of EDLC manufactured by Honda R&D. 


4.2. Computer simulation of electrochemical impedance 
with TLM 


4.2.1. Influences of various parameters on the simulated 
results of electrochemical impedance 

The results of electrochemical impedance calculated by Eq. 
(15) are shown in Fig. 10. All impedance spectra described the 
locus of lumped and distributed constant-types in low and high 
frequency ranges, respectively. Fig. 10(a) shows the simulated 
results of electrochemical impedance with variations of depths 
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Fig. 9. Experimental result of electrochemical impedance of EDLC, frequency 
range: 1 kHz—10 mHz, charging voltage: 2.0 V, and measurement temperature: 
278 K. 


of macro, meso, and micro pores. Depths of macro, meso, and 
micro pores are defined as Xma, Xme, and Xmi, respectively. 
The transition frequency shifts to low frequency side with the 
increases Of Xma, Xme, and Xmi. This result indicates that the 
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Fig. 10. Dependence of various parameters on the simulated results of electrochemical impedance of pore electrode with repeated structure (a) dependence of 
pore depth. The values of pore depth are indicated in Fig. 9(a). Dma: 50nm, Dye: 10 nm, Dmi: 2 nm, Ci 2.0 x 1074 F/cm?, p: 5.0 x 108 2cm, (O) Xma: 75 nm, 


Xme: 15 nm, Xmi: 3 nm, (A) Xma: 37.5 nm, Xme: 7.5 nm, Xmi: 1.5 nm, ( 


) Xma: 18.75 nm, Xme: 3.75 nm, Xmi: 0.75 nm and (b) dependence of diameter of a pore. 


The values of pore depth are indicated in (b), Xma: 75 nm, Xme: 15 nm, Xmi: 3 nm, cy :2.0 x 1074 F/cm?, p: 5.0 x 108 2cm, (O) Dma: 50 nm, Dme: 10 nm, Dmi: 


2nm, (A) Dma: 100nm, Dme: 20nm, Dmi: 4 nm, (l 


) Dma: 200 nm, Dme: 40 nm, Dmi: 8 nm and (c) dependence of specific resistance. The values of pore depth 


are indicated in (c), Xma: 75 nm, Dma: 50 nm, Xme: 15 nm, Dme: 10 nm, Xmi: 3 nm, Dmi: 2 nm, cy :2.0 x 1074 F/cm?, (O) 5.0 x 108 2cm, (A) 2.5 x 108 Qem, 


(0) 2.5 x 107 Qem. 
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Fig. 11. Dependence of the equivalent circuit of electrode/solution inter- 
face on the simulated results of electrochemical impedance of pore electrode 
with repeated structure changed. (a) The case that equivalent circuit includes 
Faradic process and consists of R** and Cii- The value of p is written in 
Fig. 10(a), Xma: 75nm, Dma: 50nm, Xme: 15nm, Dme: 10nm, Xmi: 3 nm, 
Dmi: 2nm, C3 : 2.0 x 1074 F/cm?, @: 5.0 x 108 2cm, a: 2.5 x 108 2cm, 
and W: 2.5 x 107 Qcm. (b) The case that equivalent circuit exclude faradic 
process, Xma: 75 nm, Dma: 50 nm, Xme: 15 nm, Dme: 10 nm, Xmi: 3 nm, Dmi: 
2nm, Re* : 6.0 x 10* 2cm, Ci : 2.0 x 1074 F/cm?, O: 5.0 x 108 Q cm, A: 
2.5 x 108 Q cm, and O: 2.5 x 107 Q cm. 


current distribution becomes non-uniform since the solution 
resistance between top and bottom of pores increases with the 
increase of X. In addition, the absolute value of impedance in 
lumped constant-type range decreases with the increase of X 
because surface area of inner wall of pores increases. Fig. 10(b) 
shows the simulated results of electrochemical impedance with 
the variations of diameter of macro, meso, and micro pores. 
Diameter of macro, meso, and micro pores are defined as Dma, 
Dme, and Dmi, respectively. The transition frequency shifts 
to low frequency side with the decrease of D because the 
solution resistance in the pores increases with the decrease 
of the cross-section area. Furthermore, the absolute value of 
impedance in the lumped constant-type range increases since 
the area of the electrode/electrolyte interface decreases with the 
decrease of D. Fig. 10(c) shows simulated results of electro- 
chemical impedance with the variations of specific resistance 
p. The transition frequency shifts to low frequency side with 
the increase of p because the solution resistance in the pores 
increases. The value of imaginary part does not change depend- 
ing on the p because the solution resistance influences on the 
high frequency impedance. Fig. 11 shows simulated results of 
electrochemical impedance of different interfacial impedance 
at electrolyte/electrode interface. Open symbols (group (a)) 
denote the impedance when the interfacial impedance is com- 
pose of the only capacitance Cf, and solid symbols (group 
(b)) denote the impedance when the interfacial impedance has 
the time constant of reaction resistance R¥* and capacitance 
C3. The impedance spectra of the groups (a) and (b) agree 


in the distributed constant-type range. Contrary to this, the 
impedance spectra of groups (a) and (b) in lumped constant- 
type range describe vertical straight line and semicircle, respec- 
tively. These results mean that the interfacial impedance of elec- 
trode/electrolyte interface Z“* controlled the shape of impedance 
spectrum in the lumped constant-type range. It was found from 
the simulated results in Figs. 10 and 11 that the interfacial 
impedance and pore dimension control the impedance in the 
lumped constant-type range and that the p and pore dimen- 
sion control the impedance in the distributed constant-type 
range. 

Simulated results of electrochemical impedance with vari- 
ous number of b are shown in Fig. 12. In the case (a) that 
bı =4, b2=16, the impedance describes the locus of lumped 
and distributed constant-types in low and high frequency ranges, 
respectively and almost agree with that in Fig. 5. However, 
the real axis and transition frequency is unclear. In the case 
(b) that bı =10, b2 =100, the transition frequency cannot be 
distinguished and the angle of straight line in the high fre- 
quency range is smaller than 45°. In the case (c) that bı = 100, 
b2 = 10,000, the locus of semicircle appears in distributed con- 
stant range. The reason is that the time constant of impedance 
of micro pore Zmi is small value since the solution resistance 
from electrode surface to the micro pore (Rsol,1 + Rsol,2) is 
small. 


4.2.2. Curve-fitting of electrochemical impedance of EDLC 
On the basis of the simulated results in Section 4.2.1, the 
curve-fitting to electrochemical impedance of EDLC was carried 
out. The fitted result is shown in Fig. 13 with the experimen- 
tal results. The fitting parameters are as follows: Xma=75 nm, 
Dma=50nm for macro pore, Xme=15nm, Dme=10nm for 
meso pore, Xmi = 3 nm, Dmi =2 nm for micro pore, Ci" = 3.0 x 
1075 F/em?, p=2.5 x 10? Qem, bı =2, bo =10, and number 
of pores is 2.0 x 10!’. The resistance of solution bulk Roulk 
was assumed as 2.2 x 107? Q, and was added to the calculated 
impedance. The simulated result of electrochemical impedance 
with this method is in good agreement with the experimen- 
tal one in Fig. 13. The diameter of a pore D is generally that 
macro pore over 50nm, meso pore 50-2 nm, and micro pore 
under 2nm. Song et al. [20] carried out the curve-fitting to 
impedance spectrum of EDLC with mode radius 0.33 wm of 
pore distribution by TLM-PSD. Since the pore size of activated 
carbon has wide distribution with a few modes models of nm 
order [6,7], the electrode pore structure of EDLC was expressed 
by three kinds of pores of nm orders in the present analysis. 
The capacitance per unit area C7" was 3.0 x 1075 Fem in 
this simulation, and this value was similar to reported ones 
3.0 x 1075 F em~? by Srinivasan et al. [11], 3.0 x 1075 F cm~? 
by Candy et al. [10] and 2.5 x 1075 F cm7? by Song et al. [20]. 
It was assumed that the specific resistance of electrolyte was 
66.7 Qcm in this simulation. This value is quite huge compar- 
ing with the specific resistance of general organic electrolyte 
because the ion migration is restricted in the pore whose size is 
1—10 nm and the electric double layer also influences on the ion 
migration in particular. In this simulation, the specific capac- 
itance p in all pores was constant value, and more accurate 
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Fig. 12. Simulated results of electrochemical impedance of pore electrode with branch structure. The numbers of branch b are (a) 4, (b) 10 and (c) 100. Xma: 75 nm, 
Dma: 50 nm, Xme: 15 nm, Dme: 10 nm, Xmi: 3 nm, Dmi: 2 nm, C7 : 2.0 x 1074 F/cm?, p: 5.0 x 108 2 cm, frequency range: (a-c) 10 mHz-10 kHz. 


analysis can be carried out to assume the different p for each 
pore. 

The merit of the present analysis is that the impedance spec- 
tra for macro, meso, micro pores can be separated in course 
of the curve-fitting to impedance spectrum. Furthermore, resis- 
tance Ry and capacitance Cy for each pore are determined by 
following formulas: 


Zm = Z4 jz% (23) 

Rez (24) 

Go (25) 
M= oZy 


In Eqs. (20)-(22), the subscript M represents the kind of pore; 
M=ma, me, and mi means macro, meso, and micro pores, 
respectively. The Rm, Cm calculated by Eqs. (20)-(22) are 
shown in Fig. 14. In Fig. 14(a), the resistance is constant value at 
low and high frequency range, and decreases with the increase 
of frequency in the transition frequency range from 0.1 to 1 Hz. 


This transition frequency in Fig. 14(a) almost agreed with that 
in Fig. 13. The order of resistance was Rmi > Rme > Rma in high 
frequency range because Rma corresponds to the bulk solution 
resistance Rpyk and Rme and Rmi includes the solution resis- 
tance in the pores in high frequency range. Resistance in low 
frequency range is larger than that in high frequency range 
in Fig. 14(a). For example, Rma involves not only Rpux but 
also the solution resistance in macro pore in low frequency 
range because the current flows to all inner wall of macro 
pore for charge and discharge. The frequency-dependent capac- 
itance for each pore is shown in Fig. 14(b). The capacitance 
has maximum value in low frequency range, and decreases 
with the increase of the frequency. The order of capacitance 
is Cma > Cme > Cmi in low frequency range. The ratio of low 
frequency capacitance of macro, meso, and micro pores is 
27.05:2.49:1, and the ratio of surface area is 31.25:2.50:1, 
which is calculated from fitting parameters for each pore. There- 
fore, it can be considered that the low frequency capacitance 
is in proportion to the surface area with regardless the kind 
of pore. 
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Fig. 13. Experimental and fitted results of electrochemical impedance of EDLC. 
The parameters for curve fitting are as follows. Xma: 75nm, Dma: 50nm, 
Xme: 15nm, Dine. 10nm, Xmi: 3nm, Dmi: 2nm, C7" : 3.0 x 1075 F/cm?, p: 
2.5 x 108 cm, bme: 10, bmi: 50, and number of pore: 2.0 x 10". 
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Fig. 14. The result of resistance, capacitance, and time constant calculated from 
Zma, Zme, and Zmi. (a) The results of relationship resistance and frequency: 
(O) Rma, (A) Rme, and (O) Rmi; (b) the results of relationship capacitance and 
frequency: (O) Cma, (A) Cme, and ( ) Cmi- 


5. Conclusions 


A new electrode model with TLM was proposed for EDLC. 
The impedance of EDLC described locus of blocking electrode 
in low frequency range and the curve at an angle of almost 45° 
to real axis in high frequency range on the Nyquist plane since 
current distribution becomes non-uniform in high frequency 
range by porous structure of electrode surface. Computer simu- 
lation of electrochemical impedance with the present electrode 
model was carried out. The contribution of five parameters to 
impedance spectrum was discussed, i.e., depth of pore, diam- 
eter of pore, specific resistance, the equivalent circuit of elec- 
trode/solution interface, and branch number. The branch number 
b influenced on the shape of impedance spectrum in distributed 
constant-type range significantly. On other hand, the interfacial 
impedance at electrode/solution interface controlled the shape of 
the impedance spectrum in lumped constant-type range. Depth 
of pore X, diameter of pore D, and specific capacitance p influ- 
ence the transition frequency. In addition, X and D influence 
absolute value of impedance in lumped constant-type range. On 
the other hand, p does not influence absolute value of impedance 
in lumped constant-type range. The resistance and capacitance 
for each pore was determined as a function of frequency, 
and the efficiency of charge and discharge in each pore was 
discussed. 
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